UNCLASSIFIED 


AD  NUMBER 


AD855766 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Jun  1969. 
Other  requests  shall  be  referred  to  the 
Air  Force  Office  of  Scientific  Research, 
Attn:  SRGO-1,  Bolling  AFB,  Washington,  DC. 


AUTHORITY 


AFOSR  ltr,  12  Nov  1971 


THIS  PAGE  IS  UNCLASSIFIED 


AD855766 


AFOSR  Final  Scientific  Report 
AFOSR-69-1399TR 


TECHNICAL  REPORT  NO.  724 


IGNITION  AND  COMBUSTION  OF  LIQUID 
HYDROCARBON  FUELS  IN  SUPERSONIC  FLOW 


Simon  Slutsky 
Forman  Williams 
Jose  Tamagno 


June  1969 


Contract  F44620-68-C-0045 
Project  9711-01 


"Research  sponsored  by  Air  Force  Office  of 
Scientific  Research,  Office  of  Aerospace 
Research,  United  States  Air  Force."  Research 
monitored  under  the  technical  supervision  of 
Dr.  Bernard  T.  Wolf son,  AFOSR.  (  $  t’PP) 


j  This  document  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 


be  made  only  with  prior  approval  of  AFOSR  (SRG^l) 

10 


u  DC 

np 

AUG  1 1989 

UlibcJLaU  U  E* 
C 


AFOSR  Final  Scientific  Report 
AFOSR-69-1399TR 


9 


TECHNICAL  REPORT.  SO , .  724  . 


IGNITION  AND  COMBUSTION  OF  LIQUID 
HYDROCARBON  FUELS  TN  SUPERSONIC  FLOW 


Simon  Slutsky 

Forman  Williams 
Jose  Tamagno 

June  1969 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Office  of  Aerospace  Research 
United  States  Air  Force 
Washington,  D.  C. 


Prepared  by 


General  Applied  Science  Laboratories,  Inc. 
A  Subsidiary  of  The  Marguardt  Corporation 
Merrick  and  Stewart  Avenues 
Westbury,  L.  I.,  New  York 


President 


S.  This  document  is  subject  to  speoial  export  controls  and  eaoh 
transmittal  to  foreign  ,-over-ments  or  foreign  nationals  may  bs 
Bade  only  Kith  prior  approval  of  AFOSR  (SRSDf{ 


ABSTRACT 


An  experimental  study  was  carried  out  to  investigate  the 
combustion  of  a  liquid  hydrocarbon  (hexane)  injected  into  a 
supersonic  stream.  Injection  modes  included  wall  and  central 
injection  configurations.  Extremely  rapid  and  intense  heat 
release  phenomena  were  obtained.  It  was  concluded  that  this 
mode  of  combustion  could  be  made  to  compete  with  that  of  pure 
gas  phase  combustion  for  many  technological  applications. 

A  preliminary  theoretical  model  of  the  combustion  mech¬ 
anism  was  constructed  which  accounts  for  many  of  the  observed 
phenomena,  and  which  explains  the  very  rapid  combustion 
characteristics . 
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PART  I 


Ry  Simon  SlufcsVy  and 
Jose  Tamagno 


I.  INTRODUCTION 

The  principal  effort  of  this  research  to  date  has  gone 
into  the  study  of  techniques  for  understanding,  predicting, 
and  then  controlling  the  processes  of  gas  phase  combustion  for 
a  supersonic  or  high  subsonic  flow.  In  the  course  of  such 
studies  it  was  observed  how  crucial  the  various  gas  transport 
processes  are  in  the  analysis  of  the  combustion  mechanism. 

Thus  it  was  concluded  that  the  rate  coefficients  of  the  various 
elementary  gas  phase  reactions  were  coupled  with  the  generally 
turbulent  diffusion,  heat  conduction  and  viscosity  character¬ 
istics  of  the  flow,  and  that  together,  they  determined  whether 
or  not  the  mixture  would  react  chemically,  and  the  distances 
characterizing  ignition  and  complete  combustion. 

It  is  the  purpose  of  the  present  study  to  extend  the  philo¬ 
sophy  of  the  previous  efforts  to  the  combustion  processes 
characterizing  high  speed  liquid  fuel-air  systems  and  ultimately 
to  more  general  multiphase  fuel  oxidizer  systems  (e.g.,  slurries, 


gels,  etc. ) . 


Such  combustion  processes  could  be  of  considerable  signif¬ 
icance  in  many  advanced  hypersonic  airbreathinq  vehicle 
applications,  especially  where  the  structural  cooling  function 
of  the  liquid  fuel  is  not  essential  for  mission  feasibility  and 
the  system  is  volume  limited.  The  advantages  of  compact  light¬ 
weight,  low  pressure  fuel  storage  facilities,  unemcumbered  by 
heat  exchanger  or  vaporizer  weight  penalties  makes  the  subject 
of  direct  liquid  fuel  injection  for  high  performance  systems  one 
of  great  interest  at  flight  Mach  numbers  between  3  and  8  and  at 
altitudes  between  60,000  and  150,000  feet. 

Such  a  study  must  involve  the  formulation  of  analytical 
models  to  describe  the  interaction  between  the  liquid  fuel 
droplets  and  the  surrounding  supersonic  airstream.  The  analytical 
representation  of  the  turbulent  mixing  and  finite  rate  combustion 
in  high  speed  flows  is,  for  the  homogeneous  gas  phase  case,  quite 
far  advanced  (Reference  1) .  Substantial  progress  has  been  made 
as  well  in  the  analytical  description  of  the  turbulent  mixing 
and  combustion  of  gas  flows  containing  cryogenic  evaporating 
liquid  components  (Reference  2) .  The  analytical  goal  of  the 
present  study  is  to  extend  this  work  by  adapting  the  previous 
experience  to  the  two-phase  (gas-liquid)  system  of  interest. 

This  involves  the  extension  and  adaptation  of  existing  mixing 


2 


and  combustion  analyses  to  the  case  of  liquid  fuel  injection, 
and  the  incorporation  of  qeneraxxZod  luuuoxti  ox  uxfxvivxcii 

and  combustion.  A  second  task  of  the  study  is  to  investigate 
experimentally  the  combustion  process  associated  with  liquid 
fuel  injection  into  a  supersonic  airstream. 

II.  EXPERIMENTAL  PROGRAM 

The  first  phase  of  the  experimental  program  was  a  preliminary 
attempt  to  investigate  the  feasibility  of  achieving  combustion 
of  a  cold  liquid  fuel  injected  into  a  supersonic  stream.  In  par¬ 
ticular,  we  started  with  an  experimental  configuration.  Figure  1, 
in  which  liquid  hexane  was  injected  normally  from  the  wall  into 
a  square  test  section.  The  fuel  flow  distribution  was  made  uni¬ 
form  across  the  duct,  insofar  ae  possible,  so  that  instruments 
and  photographs  could  be  interpreted  two-dimensionally. 

Figure  1  is  a  cross-section  of  the  test  section  including 
the  viewing  windows  and  the  wall  injectors,  and  Figure  2  is  a 
schematic  longitudinal  section  showing  the  Mach  2  nozzle,  the 
test  section,  and  an  extension  section  in  which  pressure  taps 
were  installed  at  one  inch  spacing. 

Combustion  was  achieved  in  this  apparatus,  with  remarkably 
little  difficulty.  In  fact,  the  heat  release  was  so  rapid, 
and  local  temperatures  so  high  that  the  main  problem  was  the 
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FIGURE  1  -  SUPERSONIC  COMBUSTION  OF  LIQUID  HEXANE  -  FLUSH  MOUNTED  NALL  INJ 


OF  EXPERT 


burnout  of  midstream  instrumentation.  This  loss  of  instrument¬ 
ation  sc*  tihst  tihs  dscisicn  v' 2 2  ir.sdc  tc  dclsy 

this  type,  of  measurement. 

Figure  3  is  a  photograph  of  a  flame  resulting  from  a  very 
low  rate  of  fuel  mass  flow  injected  along  one  wall,  correspond¬ 
ing  to  4  percent  (to  *  .04)  of  the  stoichiometric  fuel  air  ratio. 
Figure  4,  by  contrast,  is  an  example  of  very  intense  combustion 
in  a  configuration  with  high  rates  of  fuel  flow  (to  *  .442) 
injected  along  both  walls.  Figure  5  shows  the  pressure  distri¬ 
bution  measured  along  one  wall.  The  lowest  curve  was  made  without 
any  fuel  injection,  whereas  the  middle  and  upper  curve  are  for 
equivalence  ratios  of  .183  and  .211  respectively.  One-dimensional 
calculations  of  fuel  air  mixtures  corresponding  to  these  condi¬ 
tions  show  that  Test  26  would  not  choke,  but  that  Test  25  should. 
It  is, therefore , very  interesting  to  note  that  in  Test  25  there 
seems  to  be  only  a  small  amount  of  upstream  influence  of  fuel 
injection  and  combustion  on  the  wall  pressure  distribution  (as 
expected) .  whereas  the  pressure  distribution  of  Test  25  shows 
that  a  large  upstream  region  is  affected.  The  one-dimensional 
flow  calculation  for  Teut  26  conditions  also  indicate  an  equili¬ 
brium  p/pt^  which,  when  corrected  by  adding  the  no  burning 
(Test  23)  rise  of  about  0.04,  caused  by  boundary  layer  growth 


P  ■  115  psia 
o  *> 

T  *  3200  R 

o 

Injector:  Flush  mounted  to  bottom  wall 
Fuel  Mass  Flow:  .03  lb/sec 
0  =  .040 

Air  Mach  Number:  2.0 


FIGURE  3  -  SUPERSONIC  COMBUSTION  OF  LIQUID  INJECTED  HEXANE 


TEST  23  TEST  25  TEST  26 

To r.a  1  Pressure:  125  psia  Air  Total  Pressure:  123  psia  Air  Total  Pressure;  125  psia 

Total  Temperature:  3000°R  Air  Total  Temperature:  3100°R  Air  Total  Temperature :  3 L20°R 


PRESSURE  RISE  INDUCED 


to  station  60",  results  in  a  theoretical  pressure  corresponding 
to  100%  combustion.  This  theoretical  value  of  p/p*.  ■  .32  is 

-remarkably  close  to- that  being  approached  by  the  distribution 
of  Test  26. 

The  large  pressure  rise  from  station  26  (inches)  to  about 
38  for  Test  26  suggests  a  strong  pressure  interaction  mechanism 
between  the  region  of  intense  heat  release  and  the  surrounding 
main  air  flow.  This  conception  is  discussed  below  in  connection 
with  the  central  injection  system.  The  wall  pressure  rise  is 
thus  seen  to  be  quite  rapid,  but  so  smooth  t\iat  the  occurrence 
of  extensive  strong  shock  structures  can  be  considered  unlikely. 

Tests  were  then  carried  out  with  this  configuration  to 
determine  the  effect  of  air  temperature  on  ignition  character¬ 
istics.  Combustion  was  obtained  with  values  of  total  temperature 
as  low  as  2700°R.  However,  when  the  total  temperature  was  fur¬ 
ther  decreased  to  2500°R,  no  flame  generation  was  observed. 
Instead,  a  photograph  (Figure  6  )  of  the  fuel  injection  region 

showed  an  initial  region  of  luminosity  which  had  a  sharp  front 
at  the  line  of  injection  and  which  faded  rapidly  as  the  flow 
moved  downstream,  and  disappeared  in  a  distance  of  about 
1-1/2  inches.  It  is  believed  that  this  photograph  represents 
flame  ignition  and  quenching  and  will  be  discussed  below. 
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In  the  next  phase  of  this  exploratory  sequence,  a  test 
series  was  carried  out  with  a  midstream  injector.  Longitudinal 
and  cross-sectional  views  of  this  configuration  are  shown  in 
Figures  7  and  8.  A  piloted  midstream  injector  was  also  designed 
for  low  temperature  operation  (Figure  9),  but  use  was  postponed 
until  later. 

Very  rapid  combustion  was  achieved  in  the  system  with  air 
total  temperature  in  the  neighborhood  of  3000°R.  In  fact,  heat 
transfer  rates  in  the  neighborhood  of  the  injector  were  so  great 
that  several  injector  rigs  were  burnt  out,  and  configurations 
with  cooling  ultimately  had  to  be  built  and  installed.  A  typical 
photograph  of  the  system  is  shown  in  Figures  10a  and  10b. 

Several  resulting  wall  pressure  distributions  are  shown  in 
Figure  11.  A  no-burning  case,  Test  40,  was  run  in  order  to 
obtain  an  estimate  of  wall  boundary  layer  and  injector  drag 
effects  on  the  pressure  distribution.  The  (spatially)  oscilla¬ 
tory  character  of  the  pressure  distribution  is  due  to  the 
presence  of  multiple  reflections  of  the  induced  pressure  wave 
field.  These  wave  structures,  although  quite  well  defined,  are 
clearly  of  low  intensity  and  may  be  regarded  as  isentropic 
compressions . 
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WOT  REPRODUCIBLE 


FIGURE  6  -  IGNITION  WITH  SUBSEQUENT  FLAME  QUENCHING 
OF  LIQUID  HEXANE  IN  SUPERSONIC  AIR  STREAM 
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Test  Sectioi 


INJECTOR  WITH  ETHYLENE-OXYGEN  PILOT 
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Pq  -  120  paia 

T  -  3000°R 

o 

Injector:  Single  Wedge  Mid-Stream 
Fuel  Mass  Flow:  .106  lb/sec 
(Q  =  .  13 


FIGURE  10a  -  SUPERSONIC  COMBUSTION  OF  LIQUID 
INJECTED  HEXANE 
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Air  Total  Pressure:  134  psia 
Air  Total  Temperature:  3000°R 
Fuel  Mass  Flow:  .071  lb/sec 
Equivalence  Ratio:  .088 

Injector  Type:  Mid-Stream  Axial  Injector 


FIGURE  10b  -  SUPERSONIC  COMBUSTION  OF  LIQUID  INJECTED  HEXANE 


Tests  42  and  34  are  very  interesting.  These  flows  are 
supersonic  throughout,  the  pressure  never  reaching  the  choking 
ratio  which  is  at  p/p  *  .5  (as  discussed  later).  An  equil- 
ibrium  chemistry  one-dimensional  stream  tube  calculation  was 
carried  out  assuming  adiabatic  walls  to  determine  the  equil¬ 
ibrium  pressure  which  the  physical  flow  should  ultimately 
approach.  This  was  done,  neglecting  the  pressure  effect  of  the 
boundary  layer  growth  and  of  the  injector  drag.  When  the  latter 
"tare"  pressure  rise  is  superimposed  upon  the  "ideal"  combustion 
pressure  rise  corresponding  to  the  air  and  fuel  flow  conditions 
of  Test  42,  it  is  found  to  agree  almost  identically  with  the 
pressure  actually  experienced.  In  addition,  a  boundary  layer 
calculation  was  carried  out  and  was  used  to  calculate  the  effect¬ 
ive  area  change  of  the  duct.  This  area  change  was  then  used  to 
calculate  the  associated  rise.  This  rise  was  found  to  be  quite 
consistent  with  that  observed  in  Test  34. 

These  are  encouraging  and  exciting  results  since  they  imply 
that  essentially  complete  combustion  of  the  liquid  hexane  fuel 
was  achieved  in  about  two  feet  of  duct  at  velocities  of  about 
3000  fps.  This  kind  of  effective  kinetic  response  is  comparable 
with  the  best  that  is  achievable  with  high  performance  gaseous 
fuels  and  is  comparable  even  to  hydrogen. 
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It  is  still  more  remarkable  that  the  static  temperature 
of  the  airstream  in  the  test  section  is  quite  low,i.e.,  about 
ISS0°R  when  ff  -  ~360<5°R  and  H  -  2.0.  Hydrogen  introduced  by 
means  of  a  parallel  flow  injector  could  not  be  ignited  under 
these  conditions  without  a  pilot  or  a  shock  inducing  flame- 
holder.  Gaseous  hexane  and  the  JP  fuels  dot  of  course,  ignite 
at  this  and  lower  temperatures,  but  ignition  delay  times  of  the 
gas  phase  hydrocarbc  is  are  normally  considered  excessively  long 
for  flight  engine  application  unless  accelerated  by  means  of 
pilot  flames. 

Clearly  the  mechanism  of  ignition  and  combustion  is  not 
simply  that  of  evaporation  from  a  droplet  surface  into  an 
oxidizing  atmosphere  followed  by  ignition.  A  model  for 
description  of  the  observed  phenomena  is  discussed  in  the  next 
section. 

Before  we  proceed  to  that  section,  however,  it  is  of 
interest  to  understand  the  mechanism  of  the  static  pressure  rise 
in  Tests  42  and  34  and  of  the  subsequent  static  pressure  falloff. 
This  mechanism  can  be  clarified  by  means  of  the  two-dimensional 
flow  field  calculation  which  had  been  carried  out  assuming  that 
hydrogen  gas,  (fuel  air  ratio  to  *  0.04)  was  injected  into  a 
Mach  2  airstream  and  was  ignited  by  a  pilot  flame.  The  heat 


20 


release  rate  from  the  hydrogen  combustion  was  first  computed 
using  a  mixing  plus  finite  rate  chemistry  program  at  constant 
pressure  (Reference  3  )  to  calculate  the  shape  of  thei  equiva- 
lent  displacement  thickness,  S*,  of  the  body.  The  concepts  of 
equivalence  between  heat  release  and  volume  displacement  as 
developed  in  Reference  4  »  was  then  used  to  calculate  a 
modified  supersonic  flow  field  (Figure  12) ,  and  the  correspond* 
ing  wall  pressure  distribution  (Figure  13) .  A  pressure  cor¬ 
rection  to  the  6*  distribution  was  then  made  and  a  correspond¬ 
ingly  corrected  flow  field  computed.  It  should  be  noted  that 
the  simplest  formulation  of  the  equivalent  volume  concept 
assumes  that  the  flow  field  affected  is  that  of  an  infinite 
free  field,  and  that  any  wall  and  viscous  effects  require 
appropriate  modification.  However,  analysis  of  the  initial 
region,  including  the  neighborhood  of  the  pressure  maximum,  is 
appropriate  within  the  framework  of  this  uncoupled  model. 

The  resulting  figures  show  the  rapid  pressure  rise  due  to 
initial  heat  release,  and  the  subsequent  pressure  drop  due  to 
re-expansion  of  the  flow  as  combustion  approaches  completion. 
Thus,  the  characteristic  dependence  of  dfi*/dx  with  ~  results 
in  maximum  pressure  rise  in  the  near  region  of  the  injector. 
Incidentally,  the  small  extent  of  the  equivalent  displacement 
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volume  should  not  be  confused  with  the  extent  of  the  mixing  region. 


which  In  large.  F  u*  the  more ,  the  two-dimensional  aspect  of  the 
flow -shows  that  the  maximum  pressure  ratio  p/pt^  oi  *3®  for  3*t 

indicates  supersonic  flow  (it  would  be  sonic  at  *  .35  for  fully 

*  t  CD 

mixed  one-dimensional  flow  with  heat  release) .  Since  the  mixing 
boundary  has  not  yet  reached  the  wall,  the  stagnation  pressure  at 
the  wall  locations  in  the  neighborhood  of  maximum  pressure  station 
is  still  essentially  that  of  the  reservoir.  Consequently,  the 
wall  flow  there  is  supersonic  at  pressure  ratios  below  about  .5. 


III.  QUALITATIVE  DISCUSSION  OF  COMBUSTOR  MECHANISMS 

It  is  believed  that  the  following  sequence  of  events,  as 
illustrated  in  Figure  14,  occurs  in  the  combustor  which  explains 
qualitatively  the  observed  ignition  and  combustion  phenomena. 

These  are  initiated  with  the  injection  of  the  liquid  fuel  streams, 
normally  to  the  airstream,  whether  from  the  wall  or  central  in¬ 
jector  arrangement.  This  stream  very  quickly  becomes  unstable 
and  breaks  into  droplets.  At  the  same  time  the  flow  must  de¬ 
celerate  to  the  stagnation  conditions  on  the  droplet  surface,  so 
that  shock  fronts  must  be  generated. 

The  relative  velocity  between  the  particle  and  airstream 
is  of  the  order  of  magnitude  of  the  injection  velocity  divided 
by  the  sine  of  the  angle  between  the  particle  direction  and  the 
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free  stream  direction.  With  this  high  relative  velocity  the 


rate  of  particle  breakup  is  extremely  high. 


At  the  same  time. 


the  minimum  time  spent  by  the  partial#  in  a.  region,  where,  the. . . 

absolute  velocity  is  equal  to  or  less  than  twice  the  initial 
velocity  estimated  from  Figure  15  is  estimated  to  be  of  the 
order  of  D/v^,  where  D  is  the  measured  penetration  distance  to 
the  point  on  the  path  where  the  particles  are  moving  at  30°  to 
the  initial  free  stream  and  v^  is  the  initial  injection  velocity. 
From  Figure  15, which  is  a  photograph  of  hexane  injection  into  a 
Mach  2  stream  under  conditions  of  the  current  tests,  we 
estimate  that  this  low  speed  residence  time  is  of  the  order  of 
one  millisecond. 

A  period  of  one  millisecond  is  quite  adequate  for  ignition 
of  a  hexane  atmosphere  provided  the  static  temperature  of  the 
gas  phase  fuel  air  mixture  is  sufficiently  high.  For  the  con¬ 
figuration  of  Figure  14  the  static  temperature  experienced  by  the 
vapor  layer  around  the  droplet  will  be  somewhat  (but  not  very  much) 
less  than  the  full  total  temperature  of  the  incoming  flow. 

In  view  of  the  foregoing,  it  is  very  important  to  study  the 
system  response  as  the  stagnation  temperature  of  the  incoming  flow 
is  decreased.  In  this  way  it  was  found  that  the  lower  stagnation 
temperature  limit  for  ignition  and  combustion  was  about  2600°  R. 


P  =  115  psia 
o 

M  =2.05 

O 

Fuel  Flow  Mass:  0.226  lb/sec 

FIGURE  15-  INJECTION  OF  LIQUID  HEXANE  USING  MID-STREAM 

AXIAL  INJECTOR 
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In  terras  of  the  preceding  hypothesis,  it.  may  be  concluded  that 
below  that  critical  temperature  the  combustible  mixture  dees  not 

remain  long  enough -within  the  region -of  high -relative  velocity— - 

and  high  static  temperature  to  experience  the  necessary  ignition 
delay  time  required  for  full  combustion.  The  photograph  of  the 
injector  during  Test  30  (Figure  6)  is  particularly  interesting 
in  this  connection.  This  shows  the  wall  injector  during  injection, 
with  the  stagnation  temperature  of  the  incoming  airstream  at  2500° R. 
The  line  of  initial  luminosity  followed  by  fading  is  interpreted  as 
ignition  in  the  region  of  highest  static  temperature  on  the  droplets 
followed  by  quenching  as  the  particles  accelerate  and  the  maximum 
surface  static  temperature  drops. 

It  should  follow  from  this  interpretation  that  optimum  com¬ 
bustor  design  will  not  be  obtained  by  injection  of  the  smallest 
possible  droplet  size  or  by  injection  of  vapor  phase.  The  contrary 
should  be  true  since  the  residence  time  at  static  temperatures  of 
the  order  of  free  stream  stagnation  will  be  very  short  for  the 
latter.  Furthermore,  in  the  case  of  droplets  there  is  always  a 
region  in  the  film  about  each  droplet  having  a  fuel-air  mixture 
ratio  with  optimum  combustibility  (minimum  ignition  delay)  which 
is  in  the  neighborhood  of  the  stagnation  point  of  each  droplet. 
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For  a  corresponding  gas  jet  this  conjunction  of  high  static 
temperature  and  optimum  mixture  ratio  can  only  occur  close  to  the 

injection  point.  Furthermore!  the  spatial  distribution  of  gas 

'  <•  ’  ------- . 

fuel  at  favorable  conditions  for  ignition  is  more  limited  than 
for  the  liquid.  Thus,  the  liquid  jet  breakup  and  droplet  forma- 

I 

tion  provides  a  larger  contact  area  with  the  environment  and  as 

•  l 

indicated  above  each  fragment  or  droplet  has  its  own  optimum 
ignition  source. 

Finally,  it  is  concluded  that  the  optimum  choice  of  liquid 
fuel  jet  diameter  should  be  such  as  to  result  in  residence  time 
in  a  high  static  temperature  environment  sufficient  to  exceed 
the  ignition  delay  time  for  the  most  favorable  gas  phase  mixture 
at  that  temperature. 

In  order  to  quantify  these  hypotheses  and  determine  the 
feasibility  of  the  proposed  mechanisms,  estimates  of  the  ignition 
characteristics  of  the  liquid  jets  have  been  made.  In  addition, 
some  theoretical  considerations  are  in  order  regarding  the  losses 
in  combustion  efficiency  associated  with  vaporization  and  with 
stagnation  pressure  losses.  These  considerations  are  undertaken 
in  Part  II  of  this  report. 


REFERENCES 


1.  Ferri,  A.,  Moretti,  G.  and  Slutsky,  S.,  "Mixing  Pro- 
ceasaa  in  Supersonic  Combustion,"  J.  Soc.  Xnd.  Appl. 

Math.,  V61.  13,  No.  1 ,  March  1565. 

2.  Edelman,  R.  B.  and  Rosenbaum,  H. ,  "Finite  Rate  Evapora¬ 
tion  of  a  Cryogenic  Propellant  in  a  Binary  Two-Phase 
Mixture  with  Application  to  the  Hydrogen-Air  System, " 

AIAA  J.  4,  No.  1,  Jan.  1966,  pp.  163-165. 

3.  Edelman,  R.  and  Fortune,  0.,  "A  Preliminary  Analysis 

of  Mixing  and  Combustion  in  Confined  Flows  with  Applica¬ 
tion  to  Ejector  Ramjet  Technology,"  GASL  TR-658,  May  1967. 

4.  Chu,  B.  T.  and  Tamagno,  J.,  "Deceleration  of  a  Supersonic 
Stream  to  Subsonic  Speeds  by  Heat  Addition,"  Part  IV  of 
AFOSR  68-0467  by  Tamagno,  Slutsky  and  Chu. 


31 


PART  II 


By  Forman  Williams 
I.  INJECTION  AND  ATOMIZATION 
.A.  Introduction 

For  simplicity,  attention  will  be  focused  initially 
on  normal  (transverse)  injection  of  a  single  cylindrical  liquid 
fuel  jet  into  a  uniform  supersonic  airstream.  Of  interest  are 
jet  penetration  distances,  droplet  size  distributions  produced 
by  atomisation,  temperature  and  gaseous  fuel  concentration  dis¬ 
tributions  in  the  vicinity  of  the  jet  and  fuel  residence  times 
in  the  high-enthalpy  region.  The  preliminary  theoretical  con¬ 
siderations  presented  in  this  section  will  refer  only  to  non¬ 
reacting  conditions;  the  connection  with  ignition  criteria  will 
be  made  in  a  later  section. 
b.  Ljfesiafeaga.  Rtyigy 

Numerous  studies  have  been  made  of  liquid  disintegra¬ 
tion  in  gas  streams.  Most  of  this  work  has  been  concerned  with 
low-velocity  conditions  under  which,  for  example,  a  large  spher¬ 
ical  droplet  will  become  oblate,  will  develop  into  an  extended 
cup-shaped  sheet  or  "bag"  with  its  rim  on  the  upstream  side, 
and  eventually  will  burst  into  a  large  number  of  small  droplets 
and  a  circular  ligament  from  the  rim  which  contracts  into  a 
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number  of  eomewhat  larger  droplets.  These  low-velocity  atom¬ 
ization  studies  are  not  directly  relevant  to  high-speed  condi¬ 
tions;  at  higher  speeds  the  disintegration  mechanism  is  quits 
different,  typically  involving  stripping  of  very  small  droplets 
f rom Ihe  surface  of  thei liquid.  There  have  been  fewer  studies  of 
high-velocity  disintegration.  It  is  of  interest  to  review  very 
briefly  some  of  the  literature  on  breakup  of  liquid  jets  in  air 
streams . 

Studies  on  techniques  of  fluid  injection  in  deLaval 

nozzles  for  achieving  thrust  vector  control  motivated  much 

1  2—6 

research,  purely  theoretical  as  well  as  experimental  with 
theoretical  interpretation,  on  injection  of  liquids  into  super¬ 
sonic  streams.  Since  the  objective  in  thrust  vector  control^is 
merely  to  produce  a  side  force,  atomization  and  penetrations  are 
not  directly  relevant.  Therefore,  although  these  studies  reveal 
many  characteristics  of  the  influence  of  injection  on  the  super¬ 
sonic  stream,  they  yield  very  little  information  on  the  behavior 

of  the  liquid  jet  itself.  Experimentally  observed  jet  penetra- 

4 

tion  distances  can  be  found  in  some  of  this  work,  but  informa¬ 
tion  on  other  jet  characteristics  was  not  obtained.  Data  on 
liquid  jet  penetration  distances  in  supersonic  streams  have 
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7  8 

been  obtained  in  other  studies  '  that  were  more  fundamentally 

8  9 

oriented,  and  theories  of  jet  penetration  have  been  presented.  ' 
One  of  these  studies,^  against  which  unfortunately  a  number  of 
objections  can  be  raised  on  technical  grounds,  was  in  fact 
motivated  by  the  supersonic  combustion  problem.  To  develop  an 
accurate  theory  for  even  so  gross  a  jet  parameter  as  penetration 

9 

distance  is  difficult,  because  complex  breakup  phenomena  influ¬ 
ence  penetration. 

Very  few  studies  have  considered  the  fuel  concentra¬ 
tion  and  temperature  fields  that  develop  in  the  vicinity  of  a 
liquid  jet  for  transverse  injection.  Apparently,  data  of  this 
type  are  available  only  for  subsonic  gas  streams.10 

Abundant  data  on  jet  breakup  mechanisms  and  resulting 

droplet  size  distributions  have  been  accumulated  for  liquid  jets 

11-18 

in  subsonic  gas  streams,  but  very  limited  experimental 

results  are  available  for  mean  droplet  sizes  and  droplet  size 

19-22 

distributions  produced  in  supersonic  streams.  A  few 

theories  have  been  advanced  for  liquid  jet  breakup  in  subsonic 

(References  9,  17-18,  23-25)  and  supersonic  (References  9,23-25) 

25 

gas  streams,  the  last  yielding  theoretical  expressions  for 
the  mean  droplet  size  in  the  atomized  spray.  The  theories 


differ  greatly  in  their  physical  assumptions  and  invoke  at 
least  three  distinct  physical  mechanisms  for  disintegration* 
as  discussed  below  (Section  l.C).  Thus,  differences  between 
theories  and  paucity  of  reliable  experimental  data  cause  esti¬ 
mates  of  disintegration  characteristics  of  liquid  jets  in 
supersonic  gas  streams  to  be  highly  uncertain. 

Some  ideas  concerning  liquid  jet  disintegration  may 

26-36 

be  derived  from  a  long  series  of  careful  experiments  on 

the  breakup  of  liquid  droplets  in  gas  streams.  Typically, 
droplets  of  controlled  size  are  produced  and  suspended  in 
shock  tubes  by  various  ingenious  techniques,  and  photographic 
observations  of  the  droplets  are  recorded  after  passage  of  a 

35.3c 

shock  wave.  Only  the  last  two  of  these  studies  achieved 

supersonic  flow  of  the  gas  behind  the  shock  wave  relative  to 
the  droplet,  and  even  in  these  two  cases  no  data  are  reported 
for  convective  Mach  numbers  based  on  relative  velocity  behind 
the  shock  in  excess  of  1.5.  However,  the  results  reveal  no 
startling  differences  between  breakup  mechanisms  for  super¬ 
sonic  Mach  numbers  and  breakup  mechanisms  for  higher  subsonic 
Mach  numbers;  droplet  flattening  into  a  lenticular  shape 
followed  by  some  type  of  surface  stripping  always  appears  to 
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occur.  Aside  from  the  minor  difference  arising  from  cylindrical 
versus  spherical  geometry,  jets  will  differ  from  droplets  in  the 
presence  of  a  convective  liquid  flow  transverse  to  the  gas  stream, 
in  the  probable  presence  of  a  much  higher  intensity  of  turbulence 
of  the  liquid,  and  possibly  in  the  presence  of  enhanced  vaporiza¬ 
tion  rates  for  volatile  liquids.  Nevertheless,  one  would  expect 
that  at  least  there  would  be  a  tendency  for  liquid  jets  to 
exhibit  the  phenomena  observed  for  droplets.  This  expectation 
is  demonstrated  by  the  fact  that  theoretical  descriptions  of 
disintegration  mechanisms  for  droplets  in  high-speed  flows  invoke 
the  same  physical  phenomena  that  have  been  used  in  discussing  jet 
disintegration  (for  example,  the  analysis  in  Reference  36  is  based 
on  a  viscous  shearing  disintegration  mechanism  whose  fundamental 
attributes  closely  reserible  those  of  the  analysis  in  Reference  17) . 

c.  Pii.k^ggm.J-pin  MssiiaaiJM 

It  is  of  interest  to  look  more  closely  at  disintegra¬ 
tion  mechanisms  for  liquid  jets  because  the  physical  mechanism 
of  disintegration  affects  the  functional  dependence  of  jet  break¬ 
up  time  on  liquid  and  gas  properties  and  because,  as  will  be 
argued  in  Section  II,  the  jet  breakup  time  is  directly  relevant 
for  ignition  criteria. 
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Three  different  physical  processes  that  can  lead  to 
disintegration  have  been  identified  clearly  in  the  -literature* 

One  mechanism,  for  which  an  analysis  was  first  given 
by  G.  I.  Taylor  in  a  study  of  the  dynamics  of  droplets  subjected 
to  air  flows,  involves  merely  the  theory  of  steady  viscous  flows, 
without  consideration  of  any  instability  phenomena.  This  mech¬ 
anism  is  analyzed  in  Reference  17  for  liquid  jeta  in  subsonia 
flows  and  in  Reference  36  for  droplets  in  supersonic  streams. 

The  flow  in  a  steady-state  two-phase  viscous  boundary  layer  on 
the  windward  side  of  the  liquid  is  described  by  simplified 
versions  of  the  boundary- layer  equations.  The  azimuthal  com¬ 
ponent  of  mass  flow  rate  in  the  liquid  phase,  induced  by  the 
surface  shear,  is  calculated  at  a  position  90  deg  from  the 
forward  stagnation  point,  and  it  is  postulated  that  because  of 
the  flattened  shape  of  the  liquid,  the  liquid  boundary  layer 
separates  from  the  bulk  of  the  liquid  at  this  point,  so  that 
the  calculated  mass  flow  rate  equals  the  rate  at  which  mass  is 
stripped  from  the  liquid  surface.  Obviously,  the  velocity  Of 
the  gas  relative  to  the  liquid  and  the  viscosities  of  both  the 
liquid  and  the  gas  will  appear  as  parameters  affecting  the 
breakup  time  in  these  theories;  surface  tension  will  be  irrel- 
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evant.  We  shall  identify  this  type  of  disintegration  mech¬ 
anism  as  the  steady-shear  mechanism. 

A  second  disintegration  mechanism  emerges  from  the 
well-known  fact  that  wind  blowing  over  a  liquid  surface  gen¬ 
erates  surface  waves.  The  flow  of  the  gas  relative  to  the 
surface  of  the  liquid  jet  may  therefore  be  expected  to  gen¬ 
erate  surface  waves  on  the  jet  by  the  same  type  of  instability 
phenomenon.  Droplet  formation  will  occur  when  the  ratio  of 
amplitude  to  wavelength  for  these  waves  becomes  too  large  for 
the  crests  to  remain  an  integral  part  of  the  jet.  Although 
the  theoretical  foundations  for  describing  this  type  of  wave 
generation  were  laid  long  ago  by  tt.  Jeffreys  and  others,  the 

application  to  jet  disintegration  was  carried  out  by 
23-24 

E.  Mayer.  Breakup  times  predicted  by  this  kind  of  theory 

will  depend  on  the  relative  gas  velocity  and  on  the  surface 
tension  of  the  liquid;  the  viscosity  of  the  liquid  enters  only 
secondarily,  as  a  damping  effect,  and  the  viscosity  of  the  gas 
does  not  enter  at  all.  We  shall  call  this  type  of  disintegra¬ 
tion  mechanism  the  capillary-wave  mechanism. 
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Analysis  of  the  third  type  of  disintegration  is  also 
traceable  to  Q.  I.  Taylor;^7  this  is  the  wall »known  "Taylor 
instability."  When  a  fluid  is  accelerated  either  by  gravity 
or  by  inertia  in  a  direction  perpendicular  to  an  interface 
across  which  there  is  a  Discontinuity  in  density,  waves  grow 
on  the  surface  provided  that  the  body-force  vector,  in  a 
coordinate  system  moving  with  the  undisturbed  interface,  is 
directed  from  the  more  dense  fluid  into  the  less  dense  fluid. 
The  upstream  face  of  a  liquid  jet  is  subject  to  breakup  by 
means  of  this  type  of  instability.  Adelberg9'25  was  the  first 
to  include  this  mechanism  in  an  analysis  of  jet  disintegration t 
it  has  never  been  applied  to  droplet  breakup.  The  resulting 
breakup  time  depends  on  the  relative  velocity  of  the  gas  and 
the  mass  of  the  liquid  but  not  on  either  surface  tension  or 
gas  viscosity;  liquid  viscosity  enters  only  in  a  relatively 
minor  way  by  producing  damping.  This  type  of  disintegration 
mechanism  will  be  termed  the  acceleration-wave  mechanism,  a 
nomenclature  which  is  more  precise  than  the  more  conventional 
"gravity-wave"  label. 
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The  analyses  referenced  above  are  necessarily  quite 
Imprecise.  It  is  understandable  that  phenomena  so  complexes 
finite-amplitude  wave  development  on  surfaces  of  complex  geometry 
cannot  be  analyzed  accurately.  The  results  constitute  at  best 
oider-of -magnitude  estimates  and  at  worst  dimensional  analyses. 
Moreover,  the  analyses  do  not  take  into  account  various  phenomena, 
such  as  vaporization  of  the  liquid  or ' turbulence  in  the  liquid 
jet,  which  certainly  will  modify  the  results  quantitatively  and 
some  of  Which  might  produce  jet  disintegration  in  the  absence  of 
the  proposed  mechanisms.  There  may  be  ways  to  account  approxi¬ 
mate  1/  for  some  of  these  phenomena;  for  example,  although  argu¬ 
ments  to  the  contrary  can  be  given,  it  may  be  reasonable  to  try 
to  descries  the  effects  of  liquid  turbulence  by  using  a  turbu¬ 
lent  viscosity  coefficient  in  place  of  the  molecular  viscosity 
of  the  liquid.  However,  for  the  purpose  of  obtaining  a  rough 
estimate  of  the  disintegration  time  for  use  in  an  ignition 
theory,  it  seems  best  to  avoid  imposing  further  complications 
on  existing  theories  and  instead  to  attempt  to  select  the  mech¬ 
anism  that  appears  most  likely  to  occur  under  the  conditions  of 
interest  in  svpersonic  combustion,  and  to  employ  as  simple  a 
resoresentation  as  \s  possible  for  describing  the  implications 
of  this  mechanism. 
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Accordingly  wa  obaerve  that  for  a  liquid  of  given 
dimension*/  as  the  velocity  of  tha  gas  relative  to  the  liquid 
increases/  the  acceleration-wave  mechanism  progressively  becomes 
relatively  more  important  than  either  of  the  other  two  mechanisms 
that  have  been  considered.  This  result  is  demonstrated  by 
Adelberg  '  in  comparison  with  the  capillary-wave  mechanism 
and  can  also  be  demonstrated  in  comparison  with  the  steady- 
shear  mechanism.  It  appears  that  under  the  high-velocity  con¬ 
ditions  of  interest  in  supersonic  combustion/  the  acceleration- 
wave  mechanism  will  be  dominant  for  liquid  jet  disintegration. 
Since  droplets  are  typically  less  massive  than  jets  and  can 
therefore  be  accelerated  more  easily/  it  would  follow  from  this 
observation  that  the  breakup  of  droplets  in  supersonic  streams 
should  best  be  treated  by  the  acceleration-wave  mechanism,  in 
spite  of  the  fact  that  Ranger  and  Nicholls^6  developed  a  steady- 
shear  theory  for  correlating  their  experimental  results.  The 
photographs  aho*'n  in  Reference  36  for  the  higher  Mach  numbers 
(convective  Mach  number  relative  to  droplets  -S  1.4)  appear  to 
demonstrate  clearly  the  development  of  rather  long  liquid 
"fingers,"  reminiscent  of  Taylor  instability ,  extending  up¬ 
stream  on  the  windward  face  of  the  droplet,  as  would  be  expected 
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for  an  acceleration-wave  mechanism  but  not  for  a  steady-shear 
mechanism.  Thus,  it  would  be  of  interest  to  attempt  to  corre¬ 
late  the  high-speed  results  in  Reference  36  by  means  of  an 
acceleration-wave  theory.  In  our  present  considerations  of 
supersonic  combustion,  we  shall  restrict  our  attention  to  the 
acceleration-wave  regime. 

-D.  Jet  Breakup  Time  and  Penetration  Length 
from  Acceleration-wave,  Theory 

The  growth  time  for  acceleration  waves  is  of  the 

order  of^7 

T  ■  Jt/iL  (1) 

g 

where  l  is  the  wavelength  of  the  disturbance  and  a  is  the 
acceleration.  A  force  balance  for  an  element  of  the  liquid 
jet  of  diameter  d  and  length  Ax  in  a  stream  of  air  (of 
density  p  )  with  relative  velocity  component  v  normal  to  the 
axis  of  the  jet.  can  be  written  as 

Am  a  ■  2”  pa  y2  CD  (^xd)  .  (2) 

where  Am  is  the  mass  of  the  element  of  the  liquid  jet  and 
CD  is  the  drag  coefficient  of  the  jet.  Since 

Am  ■  j  d2  Ax  p^  (3) 
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where  P,  is  the  density  of  tha  liquid,  we  find  from  Equations 

V 

(1)  and  (2)  that  -  . . . .  '  ” 


r 


g 


(4) 


If  t  is  of  the  order  of  d  and  if  CD  is  of  the  order  of  tr/2, 
then  Equation  (4)  reduces  approximately  to 

T  -  (d/v)  Jp./P.  '  (5) 

y  v 


This  characteristic  growth  time  for  acceleration  waves  can  be 
interpreted  as  the  jet  breakup  time  that  would  result  from  an 
acceleration-wave  breakup  mechanism. 

If  t ^  is  the  breakup  time  of  the  jet,  then  the  pene¬ 
tration  distance  x^  of  the  liquid,  measured  along  the  axis  of 
the  jet,  is  given  approximately  by 

x  ■  v jT  *  d  (v./v)  Jp . /p  ""1  ,  (6) 
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where  v^  is  the  injection  velocity  of  the  liquid.  This  formula 
neglects  such  phenomena  as  further  penetration  of  droplets  after 
atomisation,  the  change  in  v  between  the  injection  point  and  the 
tip  of  the  jet  arising  from  curvature  of  the  jet  axis,  and  the 
decrease  in  d  with  increasing  x  due  to  atomization.  The  equation 
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is  useful  only  for  obtaining  the  order  of  magnitude  of  xp  . 
Since  the  mass  rate  of  injection  through  a  single  orifice  is 
given  by 

A  -  J  d2  ,  (7) 

noticing  that  tr/4  is  of  the  order  of  unity  we  can  write 
Equation  (6)  in  the  alternative  form 

xp  -  m/(dv  )  (8) 

Equations  (6)  and  (8)  are  presented  here  only  for  general 

* 

interest  and  will  not  be  needed  in  the  ignition  analysis. 

A  more  accurate  theory  which  removes  many  of  the 
approximations  that  are  implicit  in  the  equations  given  here 
may  be  found  in  References  9  and  25.  The  improved  theory 
uncovers  a  number  of  additional  parameters  for  which  approxi¬ 
mations  must  be  introduced  due  to  the  complexity  of  the  phenom¬ 
ena  involved.  It  also  leads  to  considerably  more  complicated 
equations.  However ,  the  results  do  not  appear  to  differ  in 
order  of  magnitude. 

*  Recently  rather  extensive  data  on  jet  penetration  distances 
have  been  taken  [M.A.Koplin,  K.P.Horn  and  R.E.Reichenback, 
AIAA  J.  6,,  853  (1968)].  The  resulting  correlations  are  in 
good  agreement  with  the  present  theory  over  nearly  three 
orders  of  magnitude  in  gas-to-liquid  dynamic  pressure  ratio, 
in  the  acceleration-wave  regime. 
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E.  Mean  Droplet  Size  in  Acceleration-Wave  Regime 

For  analysing  the  two-phase  spray-combustion  processea 
occurring  after  atomisation,  one  must  know  the  drop-size  distri¬ 
bution  produced  during  atomization.  No  experimental  data  at  all 
on  size  distributions  are  available  in  the  acceleration-wave 
regime,  and  apparently  only  the  rudimentary  data  of  Reference  21 
are  available  for  mean  droplet  sizes  in  this  regime.  A  theoreti¬ 
cal  result  for  the  number -weighted  average  drop  diameter  0  in  the 
acceleration-wave  regime  is  derived  in  Reference  25;  it  is 

e  -  65.3  ^V*T/(Pav2)]2/3  (9) 

where  and  o.  are  the  viscosity  and  surface  tension  of  the 
liquid.  Many  questions  can  be  raised  concerning  the  analysis 
that  leads  to  Equation  (9).  It  is  not  clear  that  the  proper 
physics  of  the  process  has  been  taken  into  account;  a  cascade 
mechanism  of  multistage  droplet  breakup  may  be  involved  in  the 
real  atomization  process,  turbulence  in  the  liquid  jet  may  be  of 
importance,  or  a  gas-phase  shearing  off  of  instability-generated 
protuberances  may  be  dominant  so  that  the  viscosity  of  the  gas 
affects  ft.  Moreover,  at  this  stage  it  appears  that  balance 
arguments  comparable  in  simplicity  to  those  of  Section  I .  D 
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cannot  yield  a  numerically  reasonable  formula  for  6,  such  ae 
Equation  (9).  Therefore  for  the  purpose  of  supersonic  combus¬ 
tion  analyses  it  seems  best  at  present  to  troat  droplet  sizes 
and  sice  distributions  as  adjustable  parameters  that  are 
assigned  values  which  lie  in  a  reasonable  range  but  are  varied 
to  obtain  correlation  with  experimental  data  on  combustion. 

It  would  be  of  interest  to  make  experimental  observations  of 
the  atomization  process  and  of  droplet  sizes  under  theee  high¬ 
speed  conditions  as  aids  in  discovering  the  relevant  physical 
mechanisms  and  in  providing  input  data  for  combustion  analyses. 

F.  Vaporization  in  the  Vicinity  of  the  Liquid  Jet 

For  the  jet  to  participate  in  the  ignition  process, 
some  vaporization  of  the  fuel  must  occur  in  its  vicinity,  since 
the  boiling  point  of  the  fuel  typically  is  well  below  the 
"ignition  temperature."  It  is  obvious  that  some  vaporization 
will  occur  because  of  the  high  stagnation  temperature,  and  it 
is  also  obvious  that  because  of  the  complicated  three-dimensional 
flow  and  the  complex  and  unknown  atomization  processes  occurring, 
an  accurate  description  of  the  temperature  and  concentration 
fields  consequent  to  vaporization  cannot  be  obtained.  Simplified 
models  must  be  developed.  The  most  pertinent  simplified  model 
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appears  cn  physical  grounds  to  be  vaporisation  in  a  two-dimen¬ 
sional  stagnation-point  boundary  layer  ©n  a  flat  plate  of 
liquid  fuel.  For  the  purpose  of  ignition  analyses,  it  ia 
appropriate  in  describing  this  flow  field  to  neglect  chemical 
heat  release  in  the  first  approximation  and  therefore  to  employ 
a  theory  applicable  to  nonreacting  flow.  This  nonreacting 
stagnation-point  boundary-layer  problem  with  coupled  heat  and 
mass  transfer  is  relatively  simple  and  has  been  analysed  pre¬ 
viously.38  41  Results  are  available  for  tanparature  and  con¬ 
centration  profiles.  The  principal  predictions  are  that  there 
is  a  smooth  concentration  profile  in  Which  the  concentration  of 
vapor  decreases  from  a  peak  value  at  the  liquid  surface  to.  sero 
at  the  edge  of  the  boundary  layer,  and  that  the  temperature 
profile  rises  smoothly  from  the  boiling  point  at  the  liquid 
surface  to  the  stagnation  temperature  at  the  edge  of  the  bound¬ 
ary  layer.  The  significance  of  these  results  for  ignition 
theory  is  that  an  extended  region  of  gas  with  concentrations 
lying  within  the  f lammability  limits  is  exposed  to  temperatures 
thr.t  are  only  slightly  below  the  stagnation  temperature. 
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II. 

A.  Approach  to  ignition  Analysis 

It  ia  apparent  from  the  preceding  observations  that 
conditions  conducive  to  ignition  exist  in  the  vicinity  of  the 
injected  liquid  jet.  Since  an  accurate  theoretical  description 
of  this  flow  region  is  unattainable#  approximate  physical  con¬ 
cepts  are  necessary  in  ignition  analyses.  Ample  background 

material  is  available  on  the  theoretical  development  of  approx- 

42-53 

imate  ignition  concepts.  For  the  present  problem,  it 

seems  best  to  divide  the  overall  combustion  reaction  into  two 
parts— an  energetically  neutral  ignition  reaction  in  which 
radicals  are  formed  through  a  chain  mechanism#  followed  by  an 
exothermic  combustion  reaction  in  which  the  radicals  recombine. 
Then  the  ignition  reaction  will  not  affect  the  temperature  field. 
Moreover,  it  appears  reasonable  at  this  stage  to  neglect  any 
differences  between  kinetic  mechanisms  in  premixed  and  nonpre- 
mixed  systems  for  the  ignition  reaction  and  to  assume  that 
during  the  ignition  process  the  gas-phase  region  in  which  the 
fuel-air  mixture  ratio  lies  within  the  flammability  limits 
remains  sufficiently  broad  in  spatial  extent  for  diffusional 
losses  of  radicals  to  be  negligible.  The  ignition  reaction  is 
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then  unaffected  by  concentration  gradients  in  the  flow  field, 

53 

and  special  consideration  of  quenching  processes  is  not 

required.  Under  these  conditions,  one  can  ascribe  to  each 

point  in  the  combustible  gas  mixture  an  ignition  delay  time 

t(T)  which  depends  only  on  the  local  static  temperature.  If 

hydrogen  were  the  fuel,  then  one  could  calculate  t(T)  from 

available  data  on  rates  of  elementary  reactions  by  using  the 

known  kinetic  mechanism)  for  liquid  fuels  such  as  hydrocarbons 

the  elementary  rates  and  kinetic  mechanisms  are  not  well  known. 

54-55 

However,  extensive  data  is  available  and  attempts  at 

formulating  reaction  mechanisms  for  hydrocarbons  have  been 
made  and  successfully  applied  to  these  dataf6'57  Thus,  we  obtain 
the  function  r(T)  from  these  studies  and  use  this  result  along 
with  the  preceding  assumptions  to  obtain  an  ignition  criterion 
as  described  below. 

I 

B.  General  Ignition  Criterion 

An  element  of  combustible  gas,  moving  at  the  local 
mass-weighted  average  velocity,  will  be  exposed  in  the  vicinity 
of  the  injector  to  temperatures  which  change  rapidly  with  time. 
We  are  therefore  faced  with  the  problem  of  applying  ignition 
delay  data,  obtained  under  isothermal  conditions,  to  a  flow  in 
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which  the  temperature  varies.  If  it  is  assumed  that  the  over¬ 
all  ignition  reaction  can  be  approximated  as  a  one-step  zero- 
order  chemical  kinetic  process  with  reaction  rate  w(T)  ■  l/r(T), 
then  the  kinetic  equation  for  the  overall  reaction  is 

dY/dt  -  «(T)  -  1/r (T)  ,  (10) 

where  Y  is  a  reaction  progress  variable  which  goes  from  zero  at 
the  beginning  of  the  reaction  to  unity  at  its  completion.  Inte¬ 
gration  of  Equation  (10)  yields 

»^i 

j  Cdt/r  ( T )]  ■  l  (11) 

o 

as  the  equation  for  the  ignition  time  of  an  element  of  com¬ 
bustible  gas  exposed  to  a  specified  time-dependent  temperature 
T(t).  Equation  (11)  can  be  used  as  an  ignition  criterion  for 
the  present  problem  in  the  following  way. 

In  the  steady-state  flow  field  established  in  the 
vicinity  of  the  injector  and  downstream,  we  trace  the  path 
followed  by  each  element  of  combustible  gas,  from  the  time  the 
element  is  formed  until  it  leaves  the  motor.  On  the  assumption 
of  nonreacting  flow,  we  calculate  the  convective  velocity  of 
the  element  and  also  the  static  temperature  to  which  it  is 
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aub-ieeted,  nv*r  ♦■'hi*  entire  history  for  each  element?  if 
Equation  (11) is  satisfied  at  any  point  inside  the  cowtoustion 
chamber,  for  any  element  of  combustible  gas,  then  we  conclude 
that  spontaneous  ignition  occurs  in  the  combustion  chamber. 

An  ignition  criterion  thereby  is  obtained  from  Equation  (11). 

C .  Approximate  Ignition  Criterion 

Since  the  flow  field  is  too  complicated  for  accurate 
calculations  of  temperature  and  velocity  histories  to  be  com¬ 
pleted,  it  is  of  interest  to  resort  to  rough  physical  ideas 
for  making  an  approximate  application  of  Equation  (11).  It  is 
well  known  that  r(T)  is  generally  expressible  in  the  form 

T  m  A  exp  (E/RT)  ,  (12) 

where  R  is  the  universal  gas  constant,  E  is  the  (constant) 
overall  activation  energy,  and  A  is  also  a  constant,  It  there¬ 
fore  follows  from  Equation  (11)  that  regions  of  high  tempera¬ 
ture  and  regions  of  long  residence  time  are  most  effective  for 
producing  ignition;  regions  where  the  temperature  is  low  or 
the  velocity  is  high  will  produce  relatively  small  contributions 
to  the  integral  in  Equation  (11).  Thus,  combustible  gases  blown 
downstream  from  the  injected  liquid  jet  will  tend  to  be  ineffec¬ 
tive  in  producing  ignition  on  two  counts— first  because  of  their 
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reduced  residence  time  and  second  because  of  their  reduced 
temperature.  The  same  comment  applies,  to  a  lesser  extent, 
to  combustible  gases  surrounding  atomized  droplets,  because 
the  droplets  are  probably  so  small  that  they  are  accelerated 
rapidly  [Equation  (9)  yields  6  ~  20  «  d].  Thus,  we  are 

led  to  the  conclusion  (neglecting  wall  effects)  that  residence 
of  combustible  gases  on  the  upstream  side  of  the  liquid  jet 
produces  the  principal  contribution  to  the  integral  in 
Equation  (11).  A  rough  ignition  criterion  can  therefore  be 
obtained  by  applying  Equation  (11)  only  to  a  gas  element  in 
the  boundary  layer  on  the  upstream  face  of  the  liquid  jet. 

Zn  Section  I.F  it  was  argued  that  combustible  gases 
are  subjected  to  temperatures  only  slightly  below  the  stagna¬ 
tion  temperature  T°  on  the  upstream  face  of  the  liquid  jet. 

To  complete  the  application  of  Equation  (11)  to  a  fuel  element 
in  this  location,  one  need  only  estimate  its  residence  time  in 
this  region  of  high  static  temperature.  To  avoid  a  complicated 
three-dimensional  boundary- layer  calculation,  we  observe  that 
the  velocity  of  this  critical  gaseous  fuel  element  is  directed 
parallel  to  the  axis  of  the  liquid  jet  and  is  approximately  equal 
in  magnitude  to  the  velocity  of  the  liquid.  Therefore  the  resi¬ 
dence  time  of  the  critical  element  of  combustible  gas,  in  the 
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region  where  T  «  T° ,  if  f imply  thf  jot  breakup  time  derived  in 
Section  I.D.  In  view  of  Equation  (5)  ^plication  of  Equation 
(11)  therefore  yields  the  following  approximate  ignition 
criterion: 

T°  ■  t(T°)  £  (d/v)  v/p4/Pa ' 

From  the  derivation/  it  is  clear  that  Equation  (13) 
represents  an  extremely  rough  result.  Taken  literally/  the 
derivation  would  predict  that  the  (exothermic)  flame  would  be 
attached  to  the  tip  of  the  injected  liquid  jet.  Experimental 
observations  seem  to  show#  on  the  contrary/  that  the  flame 
develops  some  distance  downstream  from  the  injection  region. 

On  the  basis  of  the  present  theoretical  ideas/  this  observation 
can  be  explained  by  saying  that  the  jet  is  not  quite  sufficient 
to  cause  ignition  under  the  experimental  conditions,  but  instead 
the  continuing  ignition  reactions  in  the  vicinity  of  droplets 
atomized  from  the  tip  of  the  jet  eventually  cause  Equation  (11) 
to  become  satisfied  at  some  distance  downstream  before  the  drop¬ 
lets  leave  the  combustion  chamber.  An  alternative  explanation 
of  the  observed  behavior  would  be  that  oblique  shocks  from  the 
walls  increase  the  temperature  in  the  vicinity  of  the  droplets, 
thereby  providing  a  step  in  the  contribution  to  the  integral  in 
Equation  (11)  which  rapidly  produces  ignition. 
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A  few  aspects  of  Equation  (13)  may  contribute  to  its 
usefulness  as  a  practical  ignition  criterion  in  spite  of  the 
lack  of  rigor  in  its  derivation.  Since  T  is  somewhat  below  T° 
in  the  critical  gas -phase  ignition  zone,  the  use  of  r°  tends  to 
underestimate  the  ignition  time  and  thereby  provide  an  over- 
optimistic  ignition  criterion.  On  the  other  hand,  since  con¬ 
tinuation  of  the  ignition  reactions  in  the  vicinity  of  shed 
droplets  is  neglected  in  Equation  (13),  an  underoptimistic 
ignition  criterion  would  tend  to  be  obtained.  Thus,  cancella¬ 
tion  of  cosseting  effects  may  contribute  to  the  practical 
accuracy  of  Equation  (13). 

D.  Design  Criteria 

Equation  (13)  has  a  very  simple  interpretation  as  an 
injector  design  criterion.  It  predicts  that  only  the  injector 
orifice  diameter  influences  ignition  of  a  given  fuel  in  a  given 
air  stream.  The  equation  therefore  can  be  used  as  a  formula  for 
the  critical  orifice  diameter  that  assures  ignition.  In  design, 
first  the  orifice  diameter  can  be  chosen  according  to  this 
formula,  then  other  injector  characteristics  (e.g.,  the  injection 
rate  per  injector  element)  can  be  chosen  in  order  to  achieve 
other  desirable  aspects  of  injector  performance  (e.g.,  jet  pene¬ 
tration)  . 
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E*  Comparison. with  Experiment 

It  is  of  interest  to  compare  the  predictions  of 

Equation  (13)  with  experiment.  Using  QhSL  data;  we  find  that 

typically,  in  the  experiments  where  burning  occurred, 

vVPa'  (d/v)  a.  70  4sec.  The  predicted  ignition  delay  time 

for  the  appropriate  test  conditions  is  r°  w  10  |isec.  Hence, 

ignition  is  predicted  and  observed.  At  a  Mach  number  of  2.13, 

Equation  (13)  predicts  that  the  minimum  stagnation  temperature 

for  ignition  is  approximately  2600° R  for  the  GASL  experiments. 
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Tn  the  experiments  of  Mestre  and  Viaud,  although  the  injection 
geometry  is  considerably  more  complex,  we  can  use  the  wall-injec¬ 
tor  orifice  size  to  estimate  that  Jp./p  1  (d/v)  10  psec  and  the 

v  £1 

stagnation  temperature  to  show  that  r°  <  1  ^sec,  thereby  again 
ob-aining  agreement  between  the  prediction  of  Equation  (13)  and 
the  experimental  observation  that  ignition  occurred. 

Finally,  it  should  be  emphasized  that  if  improved 
accuracy  is  needed  or  if  other  injection  geometries  are  of  in¬ 
terest,  Equation  (11)  should  be  used  in  place  of  Equation  (13). 


III.  SPRAY  COMBUSTION 


Although  injection  and  atomization  proceasea  are  ao  complex 
that  the  relevant  conaervation  equationa  cannot  be  aolved  numer¬ 
ically,  it  is  quite  possible  that  the  spray  combustion  processes 
occurring  downstream  from  the  atomization  region  can  be  analyzed 
accurately  by  such  numerical  approaches.  The  relevant  conservation 
equations  for  steady  flows  have  been  given  in  Reference  59.  To 
include  the  gas-phase  ignition  process  in  the  analysis  may  prove 
difficult,  but  in  the  region  downstream  from  the  position  at 
which  ignition  ia  completed,  it  should  be  possible  to  treat  a 
single  heterogeneous  combustion  reaction  and  thereby  simplify  the 

calculations.  Much  theoretical  work  remains  to  be  done  in  analyz- 

57 

ing  how  the  mathematical  properties  of  the  spray  equation 
modify  the  usual  characteristic  calculations  for  supersonic  flow; 
studies  of  this  kind  should  be  completed  before  attenpts  are 
made  to  program  the  full  equations.  However,  there  appears  to  be 
nothing  in  principle  that  would  prevent  the  performance  of  calcu¬ 
lations  with  the  full  equations. 

In  view  of  the  probable  existence  of  nonuniform  fuel  distribu¬ 
tions,  it  does  not  seem  desirable  to  restrict  one's  attention  to 
qua si -one -dimens ional  forms  of  the  equations,  as  was  done  in  the 
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application*  given  in  Reference  57.  At  the  very  least,  char¬ 
acteristic  analyses  for  steady-state  two-dimensional  flows 
should  be  performed.  However,  other  approximations,  such  as 
the  assumption  of  a  monodisperse  spray  and  the  assumption  that 
the  droplets  and  gas  travel  at  the  same  velocity,  may  introduce 
useful  simplif ications  without  impairing  accuracy  severely. ^ 
Work  on  the  solution  of  steady-state  two-dimensional  conserva¬ 
tion  equations  for  sprays  must  certainly  be  initiated  if  useful 
calculations  of  combustion  efficiency  are  desired. 

A  very  rough  approximation  for  combustion  efficiency  can 
be  obtained  by  substituting  suitable  numbers  into  Equation  134) 
of  Reference  59.  For  a  velocity  of  4000  ft/sec  and  a  eorrtbustion 
chamber  length  of  4  ft,  if  the  square  of  the  droplet  diameter 

decreases  linearly  with  time,  then  for  a  burning  rate  constant 

-2  2 

of  10  cm  /sec,  an  initial  mean  droplet  diameter  of  30  p  and  a 
Rosin-Rammler  3ize  distribution  with  exponent  4,  we  find  a  com¬ 
bustion  efficiency  of  9896. 
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XV.  STAGNATION  PRESSURE  IQSSgg 


Calculation  of  loaaaa  in  atagnation  pressure  should  be 
performed  for  comparing  liquid-injection  and  gaseous-injection 
supersonic  combustion.  In  the  absence  of  such  calculations, 
preliminary  thought  leads  one  to  believe  that  the  losses  asso¬ 
ciated  with  heat  release  in  supersonic  flow  are  likely  to  be 
comparable  in  the  two  cases.  However,  the  probable  lower  rate 
of  heat  release  (after  ignition)  with  droplet  combustion  may 
provide  the  designer  with  greater  opportunities  for  changing 
motor  geometry  to  alleviate  losses.  Also,  it  is  quite  clear 
that  the  stagnation  pressure  losses  associated  with  injection 
of  a  liquid  fuel  are  much  less  than  those  associated  with 
injection  of  the  same  mass  flow  rate  of  gaseous  fuel.  There¬ 
fore,  from  the  viewpoint  of  reducing  the  stagnation  pressure 
losses,  use  of  liquid  injection  may  prove  to  be  attractive. 
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